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Abstract 

The Minimum Rate of Dissipation Principle (MRDP) affirms that, for time- 
independent boundary conditions, a thermodynamic system evolves towards a 
steady-state with the least possible dissipation. In this note, examples of diffusion 
processes of two solutes in an isothermal system are analyzed in detail. In par- 
ticular, we consider the relaxation of the system when the metric tensor (i.e. the 
Onsager matrix) is constant and when the Onsager coefficients weakly depend on 
the spatial derivatives of the concentrations. We show that, to leading order, dur- 
ing the relaxation towards a steady-state, the system traces out a geodesic in the 
space of thermodynamic configurations, in accordance with the MRDP. 



1 Introduction 



The Minimum Rate of Dissipation Principle (MRDP) describes general thermody- 
namic systems relaxing towards stationary non-equilibrium states (see refs Jl|]-||2D)- 
An interesting application of this principle is to the study of magnetically confined 
plasmas relaxing towards a steady state in the highly collisional regime (Pfirsch- 
Schliiter regime) and in the low-collisional regimes (the Banana and Plateau regimes). 
However, before accomplishing this task, we will first test the validity of this prin- 
ciple by analyzing the relaxation processes of several simpler thermodynamic sys- 
tems. A good experimental test would be to analyze relaxations of solar storms 
through satellite experiments and check whether the MRDP is really satisfied dur- 
ing the evolution. Theoretical checks can be performed by analyzing the relax- 
ations of thermodynamic systems towards a steady-state through the constitutive 
evolution equations. This note is devoted to the study of two concrete examples of 
relaxations. 

Let us now briefly describe the relaxation of a thermodynamic system to a 
steady-state. In the thermodynamic field theory (TFT) description of a thermo- 
dynamic system a homogeneous configuration corresponds to a point x in the 
thermodynamic space M. The corresponding thermodynamic forces, which are 
completely determined by the configuration x, are assembled into the vector U 
defined as 

exp x (U) =iu(t f ) =0. (1) 

Here 7^7 is a shortest path from x to the point 0, which corresponds to the equi- 
librium configuration. The exponential map exp x in Eq. (d) associates to every 
tangent vector V in TM X a point y in M which lies a distance \ V\ along the unique 
geodesic which is parallel to V at x (see for example ref. [3]). The system begins at 
x at time and reaches the state y = 7(i/) at a finite time tf. 

The thermodynamic flows J M , which are dual to the thermodynamic forces £/ M , 
and defined to be 

J M = \^u u = (g^u + U„)u v (2) 

where A M „ is an asymmetric tensor, g, the symmetric part of A, is the Riemannian 
metric on M and the skew-symmetric part is The entropy production a at the 
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point x is defined to be 

o = (g^u + f^WU" = g^UW > (3) 

where the inequality corresponds to the second law of thermodynamics. As the 
Riemannian metric is positive definite, this inequality is always satisfied. If the 
point x is not a steady-state then the system will evolve. We consider the process 
known as relaxation in which the initial velocity vanishes and the system relaxes 
to a steady-state y. 

The TFT description of relaxation mainly rests upon the following postulate: 
The manifold (M,g) is a Riemannian manifold with metric g. The Minimum Rate 
of Dissipation Principle establishes that during the process of relaxation the con- 
figuration traces out a geodesic 7 in the thermodynamic space. In Refs. [10- A the 
geodesic property was referred to as the shortest path principle. It allows one to 
write 

x = l(U), y = j(tf) = exp x (V) (4) 

where V is a tangent vector to 7 at x whose norm is the geodesic distance from 
x to y. The system begins at x at time and reaches the steady-state y at the 
time t f . The process of relaxation begins as soon as the system is released, that 
is, immediately after time U. In support of this conjecture, we have argued that 
V automatically satisfies a manifestly covariant form of the Minimum Entropy 
Production Theorem (valid when the system is near equilibrium), the Glansdorff- 
Prigogine Universal Criterion of Evolution and the Minimum Rate of Dissipation 
Principle. 

In general, the macroscopic description of thermodynamic systems gives rise 
to state variables that depend continuously on the spatial coordinates. The Brus- 
sels school results mentioned above allow us to treat spatially extended thermo- 
dynamic systems by constructing the following fiber bundle: the ordinary space 
R 3 is combined with the thermodynamic space attached to each point in _R 3 (this 
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approach can been found in refs |4]-[6]). The ordinary space R 3 is the base mani- 
fold and the thermodynamic space is the fiber. The G-P theorem, obeyed on each 
geodesic, takes the integral form 

dV 

^dv (5) 

where dv denotes the volume element and the integration is over the entire space 
fl occupied by the system. 
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In this paper the relaxation of a distributed system is illustrated by analyzing 
two concrete examples. Section |2] describes the relaxation of a chemical system 
when the metric is flat (i.e., the Onsager matrix is constant). Then in section [3] 
the linear (Onsager) transport coefficients are corrected by terms which linearly 
depend on the thermodynamic forces (i.e., the Onsager transport coefficients de- 
pends on the derivatives of the chemical concentrations). We will see that the Min- 
imum Rate of Dissipation Principle is satisfied in both of these cases. 



2 Diffusion of Two Solutes in an Isothermal System with Con- 
stant Diffusion Coefficients 



In this section we analyze in detail a simple example of relaxation: the diffusion 
of two solutes in an isothermal system when the Onsager matrix is constant. The 
evolution of the mass fractions q is described by the differential equations I^Hsl] 

= D n A a + D 12 A c 2 

at 

^ = D 12 A Cl + A c 2 (6) 
at 

where A is the Laplacian operator. Eqs © are quite general. They describe a di- 
rect response of the plasma to the presence of a thermodynamic force. A radial 
pressure gradient produces a radial particle flux (Fick's law). The corresponding 
coefficient is the radial electron (ion) radial diffusion coefficient (i.e., the loss of mat- 
ter). However, Eqs © are also valid for describing general diffusion experiments, 
which occur e.g. in isomer mixtures or in isotopic mixtures (with not too different 
molecular masses) and mass diluted systems (i.e. c\ <C 1 and c 2 ~ 1). 

In this section we will consider the Onsager regime, in which the diffusion 
coefficients are uniform. In this case Eqs © can be diagonalized and rewritten 
in terms of the solutes that diffuse independently 

£ = *A* (7) 
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where the c$ are linear combinations of the mass fractions c,. 

Let us now suppose that our system is infinite and satisfies the initial conditions 



cx{T 1 t=0) = N 1 6(r) 



c 2 (r,t = 0) = N 2 6(r) 



(8) 



where the Ni indicate the total number of particles of species i. We will now ana- 
lyze the relaxation of this system and verify the validity of the Minimum Rate of 
Dissipation Principle. It is easy to check that the solutions of Eq. with initial 
conditions Eq. © can be cast into the form 



Notice that at each moment t the mass fractions are Gaussian distributions cen- 
tered at the origin with one standard deviation <7j equal to yJIDfi. The origin fol- 
lows a geodesic in the space of mass fractions; however, we are interested in the 
space of thermodynamic forces and we will see that the forces always vanish at the 
origin in this example. 

To understand the sense in which the evolution is a geodesic in a space of ther- 
modynamic forces, it will be necessary to follow the relaxation of a point at the 
time-dependent radius 



which lies a fixed number of standard deviations from the origin. We will see 
that the thermodynamic forces at r p follow a geodesic for every fixed value of p. 
One might object that we should instead be following a point at a fixed value of r; 
however, we will see that such points begin already relaxed, then they are excited, 
and then they relax again. Thus the shortest path principal, which applies only 
to the phenomenon of relaxation, does not apply to the value of a field at a fixed 
value of r in the present example. 

The thermodynamic forces are defined to be the gradients of the mass fractions. 





(9) 




(10) 
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x da rJVi / r 2 

U [r,t) = — — = , — -TT^-expi 



In particular their radial components are 

dc\ 

~dJ = ~ 16n 3 / 2 (D 1 t) 5 / 2 C " P V~4D^ 
and so they are related by 
At r = r p , Eq. ((12]) reduces to 

In particular the ratio of the forces U 1 and U 2 is independent of t and so, for each 
value of p, U 1 (t) and U 2 (t) lie along a line in the two-dimensional space of forces. 

We are considering the Onsager regime, and so the metric is flat and the line is 
a geodesic. Thus at each r p the system traces out a geodesic from a p-dependent 
excited state to the stationary state 0, which in this case is also the equilibrium state. 
In the Onsager regime the Glansdorff-Prigogine quantity is just the rate at which 
the forces approaches equilibrium with respect to the affine parameter, which is 
maximized when the velocity vector points towards equilibrium, as it always does 
for a straight line ending at equilibrium. Thus the minimum rate of dissipation 
principle is satisfied, as it must be for any geodesic relaxation according to the 
general arguments reported in refs (111]) and (PI). 



3 Diffusion Processes of a Thermodynamic System in which the 
Metric Tensor is a Small Perturbation of the Onsager Matrix 

We shall now analyze another example: the relaxation of a system when the On- 
sager coefficients (weakly) depend on the derivatives of the chemical concentra- 
tions. We shall show that the minimum rate of dissipation principle is satisfied in 
this case as well. For simplicity we will consider a system in which the concentra- 
tions depend on a single spatial direction, with coordinate x. 
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We now reconsider the diffusion process examined in the previous section, 
adding a small perturbation to the metric tensor g^ v with respect to the Onsager 
tensor L^ u 

g^u = V + A« v c' K (14) 

where L is diagonal, the symbol ' is the derivative in the x direction and A K V = A* 
are considered to be small constant parameters. 

The relaxation processes can be analyzed by solving the following evolution 
equations for infinite systems 

— = -V-Ji Ji = 011 VCi + 012 VC2 

^ = -V-J 2 J 2 = 012VC1 + 022VC 2 . (15) 

at 

Notice that the structure of these equations is quite general. For example Eq © 
describes diffusion, as well as the thermal processes of Ref. Jlo| . In this last case, a 
radial electron (ion) temperature gradient produces a radial electron (ion) heat flux 
(Fourier's law). The corresponding coefficient is the radial electron (ion) thermal con- 
ductivity. As we are considering uni-dimensional systems, the gradient operator is 
just V = (Jj, 0, 0) and the evolution equations take the form 

ci = DicJ + A x d x c[ + A 2 c 2 c 2 + A 3 (c[c 2 " + c[d 2 ) 

c 2 = D 2 c 2 + A±d x c[ + A 5 c 2 c 2 + A§(d x c 2 " + c'[c 2 ) (16) 

where the dot and prime indicate the time and the spatial derivatives respectively 
and D k = L kk for k = (1, 2). Moreover 

A\ = ^A\ x ; A 2 = 2A\ 2 ; A 3 = A\ x + A\ 2 ; A4 = 2A\ 2 ; A 5 = 2A\ 2 ; A 6 = A\ 2 + A 22 . 

(17) 

At the zeroth order in A, we have to solve the equations 

cio = -Dic'i'o 

c 20 = D 2 c 20 (18) 

with the boundary conditions 

c 10 (r,t = 0) = JV 1( f(r) 

c 20 M = 0) =N 2 6(r). (19) 
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We find 



Ni 



( —x 



Cio = (4„ Dl t)v* exp ywrt 

N 2 



C20 



(AnD 2 ty/ 2eXP \AD 2 tJ' 



x 



(20) 



As we shall see later on, our calculations are greatly simplified if, instead of coor- 
dinates (x, t), we use coordinates (x, a) defined as 



x = aa, t = — , D 



D1D2 
D 2 -D 1 



(21) 



so that a = x/y/ADt. Eqs (|Z0j) take then the form 



ClO 



C20 



y/nx 



exp(— Di& 



/'1 



X 



<\\JD 2 N 2 / r\ 2\ ^2 



y/nx 



exp(—D 2 a z 



x 



(22) 



where 

fj, k = aN k exp(-D k a 2 ) 

Using 



D k N k 



-x 



D k = ^r, with (A; = 1,2). (23) 
u k 



2a 3 D 



& x const 



x* 



r . la 

C \t const 

a X 



(24) 



and expanding the evolution equations (|16|) in A by inserting (|22|) one finds at 
leading order 



c k = D k c k - 8AI 



x N\N v D\D v oP 



(a 2 (D x + D v ) - 1) exp(-(D A + D v )a 2 ) (25) 



where summation over indexes A and v (with A and v taking the values 1 and 2) is 
implicitly understood. It is not difficult to check that solution of Eq. (|25|) is 



c k = c 0k -exp(-D k a 2 ) + e^N x N v D x D v —exp(-{D x + D v )a 2 ) (26) 
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where 



VD k N k 

cok — f= — 

V 71 " 

A x 

e x ku = ^. (27) 

D x + D v - D k 

We will now argue that, during the relaxation, the thermodynamic system 
evolves towards a steady-state along a geodesic. Recall that the thermodynamic 
forces are defined to be the spatial derivative of the mass fractions. As is easily 
checked, to first order in e, the derivatives of the affine parameter r and of thermo- 
dynamic forces, with respect to time t, behave as 

+ 0(e 2 ) 



dr 


1 


dt 




dd 


1 


dt 


Id* 



1 



(28) 



On the other hand, to first order in e, we also have 

which is consistent with the geodesic equation, confirming that the evolution of 
the system traces out a geodesic in the space of thermodynamic configurations. 
Thus the minimum rate of dissipation principle is satisfied, as it must be for any 
geodesic relaxation as was demonstrated in refs (IIO) and (||2j]). 



4 Conclusions and Perspectives 



The Minimum Rate of Dissipation Principle (MRDP) provides predictions for gen- 
eral thermodynamic systems relaxing towards stationary non-equilibrium states 
(see refs 111] and In principle, this task can be accomplished only after having 
set up a formalism where all variables are free and the structure of the fiber bun- 
dle is completely determined. In general, the macroscopic description of thermo- 
dynamic systems gives rise to state variables that depend continuously on space 
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coordinates. In this case, the thermodynamic forces possess values associated to 
each spatial point. To treat spatially extended thermodynamic systems, we have to 
construct a fiber bundle by combining the ordinary space B = R 3 with the thermo- 
dynamic space F attached to each point in R 3 . In this and in previous papers, we 
have treated non-homogenous thermodynamic systems supposing that the fiber is 
trivial i.e., the fiber bundle is simply the product space B x F. This is, however, a 
strong assumption, in fact there are a wide variety of examples where the spaces 
cannot be described as trivial fiber bundles. 

In this paper we have checked the validity of this principle by analyzing the 
relaxation processes of several thermodynamic systems, assuming that the fiber 
is trivial. In particular we analyzed the relaxation of two concrete examples: sys- 
tems where the metric tensor (the Onsager matrix) is independent of the thermo- 
dynamic forces and the case where the Onsager matrix weakly depends on the 
gradients of the concentrations. The structure of these equations is quite general: 
they can be used for describing diffusion, as well as thermal processes. A typical 
example is a radial electron (ion) temperature gradient, which produces a radial 
electron (ion) heat flux (Fourier's law). The corresponding coefficient is the radial 
electron (ion) thermal conductivity. We proved that to leading order the evolution 
of the system traces out a geodesic in the space of thermodynamic configurations. 
From ref. |lj], we then have also shown that in these cases the systems evolve to- 
wards a steady-state with the least possible dissipation (the Minimum Rate of Dis- 
sipation Principle). 
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